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Abstract 

We present a thermodynamically consistent mesoscopic model of protein adsorption 
at liquid-solid interfaces. First describing the equilibrium state under varying protein 
concentration of the solution and binding conditions, we predict a non-trivial (non- 
monotonic) dependence of the experimentally observable properties of the adsorbed 
layer (such as the surface density and surface coverage) on these parameters. We 
subsequently proceed to develop a dynamical model consistent with the equilibrium 
description, which qualitatively reproduces known experimental phenomena and 
offers a promising way of studying the exchange of the adsorbed proteins by the 
proteins of the solution. 
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1 Introduction 

Protein adsorption at liquid-solid interfaces is a fundamental problem of sev- 
eral diverse areas of biotechnology. A short list would include, biocompatibility 
of implants, blood clotting, filter fouling and protein chip technology. It is also 
a problem of general theoretical interest, especially because despite the con- 
siderable body of experimental data and empirical phenomena - due in part 
to the very precise methods available for the measurement of the adsorbed 
amount - the literature is thus far lacking a mesoscopic (particle level) de- 
scription that can account for all of the experimental observations [1]. Below 
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wc present a model that will - wc hope - not only meet the above criteria, 
but do so in a physically consistent manner. 

The basic treatise behind various particle level descriptions of protein adsorp- 
tion is that adsorbed proteins can be in several different states of different 
surface (or footprint) sizes. This is supported by experimental evidence that 
adsorbed proteins undergo surface-induced conformational change [2,3,4] char- 
acterized by a substantial growth of the surface contact area. Combined with 
evidence [5] that under conditions of high surface coverage these transitions 
are sterically hindered by neighboring proteins on the surface, the usual meso- 
scopic description is completed by assuming a hard-core repulsive interaction 
between proteins - both on the surface and in the solution. 

Although proteins adsorb to a two-dimensional (2D) substrate, we will further 
use a one-dimensional (ID) approximation. This simplification is justified, on 
one hand, by the fact that there is no thermodynamic phase transition (which 
is usually sensitive to the dimensionality) in the adsorption process, and on 
the other hand, because the main features of adsorption (such as jamming, 
slow convergence toward equilibrium, short-range spatial correlations) are in- 
dependent of the spacial dimensions. ID models already have the ingredients 
that make their behavior non-trivial and qualitatively similar to that observed 
in higher dimensions, and they are much easier to treat analytically [6]. 

Surface-protein systems of interest typically exhibit tightly bound protein 
states with binding energies upto several tens of [7]. This is the reason 
that most models to date have incorporated at least partial irreversibility to 
explain the obvious deviation of protein adsorption from simple Langmuir-like 
behavior [8]. The classical irreversible model of particle adsorption, random 
sequential adsorption (RSA) (for a recent review see [6]), in which adsorption 
is considered to be permanent and the adsorbed system of particles on the 
surface is eventually described by a jammed state with no more free space 
available for particle deposition, is inherently incapable of describing any pro- 
cess involving desorption of proteins. The second class of models [9,10,11,12,13] 
differentiate between reversibly adsorbed states of small size and irreversibly 
adsorbed larger states. These models rely on the assumption that the escape 
from the strongly bound states is practically impossible under reasonable ex- 
perimental time scales. This notion is supported by the experimental evidence 
of long relaxation times of the observed properties. We show that slow relax- 
ation of experimental observables (typically surface density) can be equally 
well described by a reversible model, and that such an approach has a non- 
trivial equilibrium state distribution, with several unprecedented predictions 
for the concentration and binding strength dependence of the surface density 
and the surface coverage, respectively. We also argue that exchange effects, 
widely observed for non-biological polymers [14], and to a limited extent for 
proteins [15] can only be treated in the context of a reversible model such as 
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Fig. 1. Adhesion potential of a protein, where each of its states i (characterized by 
a binding energy Sj and a footprint size li) is separated from the neighboring states 
by potential barriers. The corresponding transition rates are kf^ (defined in Sec. 3). 

ours. 



2 Equilibrium model 



In our model we consider a ID line of adsorption of size L in contact with a 
protein solution of constant concentration c. Each proteins being adsorbed on 
this line is considered to be in one of M different states, each characterized by 
an Ei binding energy (1 < i < M) and a monotonically increasing k linear 
footprint size (see Fig. 1). Proteins in the smallest footprint state (i = 1) can 
undergo adsorption from and desorption to the protein solution, characterized 
by a fca adsorption and a desorption rate. Adsorption and desorption in 
other states are neglected, because these states have lower energies on the 
surface and higher energies in the solution than those of the smallest footprint 

^ we can assume {li} to be monotonically increasing as for any set of {li} and 
{sj} we can obtain this by appropriately rearranging the {li}-s, because there is no 
constraint on the {ei}-s. 



3 



state, and thus, represent a much slower adsorption-desorption kinetics. Pro- 
teins interact via a hard-core repulsion, i.e., they are not allowed to overlap. 

Given a set of values for the parameters (L, c, {si}, {k}, kg., k^) the model has 
a unique equilibrium distribution of states {rij}, where rij denotes the number 
of adsorbed proteins in state i, and 

M 

T.^i^N (1) 

1=1 

is the total number of proteins adsorbed on the surface. The distribution must 
obey the 

M 

Y,nili<L (2) 
1=1 

spatial constraint. Our reversible dynamical models (defined later) will be 
explicitly shown to converge to this equilibrium distribution. 



2. 1 The conditional free energy 

The above equilibrium distribution can be obtained by determining the con- 
ditional free energy 

F({n,}) = E{{ni}) - T5({n,}) - /.iV({n,}) (3) 

for any given distribution {ui} (obeying the spatial constraint), and finding the 
distribution {n™"} which minimizes it. The energy E[{ni\) and the entropy 
S{{ni}) can be expressed as 

M 

E{{ni})=Y,^^n^ and 5(K}) = -A^b lnC(K}) (4) 

i=l 

and N is given by Eq. (1). Now we need to count the number of possible 
arrangements (configurations) C({nj}) of the proteins on the surface, and 
express the chemical potential ^ of the proteins in the solution in terms of 
their concentration and adsorption/desorption rates. Restricting ourselves to 
a ID model has the advantage that the number of configurations C({ni}) can 
be calculated in a straightforward manner. 
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2.1.1 Adsorption on a lattice 



Representing the surface as a ID lattice of period 6 (with L/S sites), the 
problem of calculating C{{ni}) reduces to a simple combinatoric exercise. 
Denoting the total occupied surface as 



M 



i=l 



and the number of unoccupied lattice sites as 

no = — (6) 



the number of different permutations of the N adsorbed proteins and the uq 
empty sites is (no + iV)!. Since proteins in the same state, as well as the empty 
sites are indistinguishable, this has to be divided by noinf^i^i' to get the 
number of different configurations: 



Considering the lattice constant S small compared to the protein footprint size 
implies that the number of empty sites is much larger than N. Consequently, 
in the above formula we can make the approximation: 

+ (8) 



no'- 



(9) 



which is shown to be exact as 5 — > 0, in the Appendix. 
Using Stirling's formula and neglecting the sub-linear terms in N: 

I L-] ^ \ 
In C({ni}) ~ \ N In — YjiMni^Ny 

The entropy can now be written as 

L - / ^ 

S{{n,}) = k^N In — — -k^Y.^, In + k^N. (10) 
" 1=1 
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Introducing the state densities 

Ui N 

i=l 



for proteins in state i and for the total adsorbed proteins, respectively, as well 
as the corresponding fractional surface coverages 

and A = = ^pjZj = -, (12) 
^ i=i i=i ^ 

we can explicitly express the extensivity of the entropy or equivalently the 
lack of dependence of the entropy density s — S/L on the system size L, by 
writing the entropy density as 



s{{Pi}) = ^^^^ = (pln(l - A) - pin 5 - f: pan A + . 



(13) 



The chemical potential can be derived from the requirement that for a sys- 
tem in equilibrium detailed balance must hold. In our case this requirement 
prescribes that the rate of adsorption to a particular site per unit time, 5ckg_, 
and the rate of desorption from the same state, k^, satisfy 

fis-ei 5ck^ 



yielding 



k 

P5 = £i + ksT In — ^ +kBT In c +kBT In 5 = p + ksT In d, 
kd 



(15) 



where we have separated p into po (which depends only on the binding con- 
ditions) and kBTlnc (which depends on the concentration). For the sake of 
brevity we will consider /cbT" = 1 throughout the paper. We can then write 
the conditional free energy density using (3), (13) and (15) as 

mPi}) = ^^^^ = Y.Piei - p(l + p) - pln(l - A) + E^lnp,, (16) 

where the In 5 terms have canceled, and hence, the free energy density is 
independent of the spatial resolution 5. 
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2.2 Equilibrium state distribution 



The free energy density (16) can be formally extremized by solving 

^^^ = £,-/. + lnp,-ln(l-A) + ^/, = 0. (17) 



This yields for the equilibrium state density distribution 

p,- = (1- A)e'^e-("^+T^'^), (18) 

which must satisfy the self-consistency criteria: 

M M 

p^Y,Pi a-nd A ^Y,hPi- (19) 

i=l 1=1 

To numerically solve the above we can proceed by defining the functions 

M 

P*(p,A) = X:P*(p,A) (20) 

i=l 



and 



M 

A*(p,A)=^Z,p*(p,A), (21) 



i=l 



where 

p*{p, A) = (1 - A) e'' e-(^*+T^'*) . (22) 

For a given set of values of the parameters (c, {si}, {/«}, k^,, kd), 

p*(p,A)-p = and A*(p,A)-A = (23) 



define two curves in the (p, A)-space, an intersection of which (at the point p = 
p"^^ and A = A''''*) gives a solution of Eq. (17). Provided that p""""^ G [0, ^^J^^ | ] 
and A*^^* G [0,1], this corresponds to an extremum of the conditional free 
energy (16), 

p^^t = (1 - A^''*) e'' e~^^*+T^''^ . (24) 
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Prom physical considerations it is obvious that there must be at least one such 
solution, which is furthermore a minimum of (16) (otherwise we would be left 
without a stable equilibrium state). Although a formal proof is so far lacking, 
numerical results indicate that there are no values of the parameters for which 
there would be more than one solution. 



2.3 Linear approximation 



If we only wish to gain information on whether, for a particular system, the 
proteins adsorbed on the surface spread out completely (energy dominated 
regime) , or whether they remain compact despite the availability of energeti- 
cally favorable but larger states (entropy dominated regime), we can proceed 
by considering a simple linear approximation for the binding energies as a 
function of the footprint size, and neglect all higher order terms: 

Si = sli. (25) 



It is immediately apparent that Eq. (24) specifies an exponential form for the 
equilibrium state density distribution: 



~(e H 

p. = (1-A)e^e =be''^^. (26) 



This result of an exponential state density distribution then, depending on the 
exponent a, corresponds either to a system - for negative a - where proteins 
remain compact compared to their energetically more favorable (assuming 
e < 0) spread out state, or - for positive a - to a system where they spread 
out and predominantly occupy their lowest energy state. Within the context 
of our approximation the crossover between these two modes of behavior can 
be analytically derived. Since for a — Q, pi ~ b the self-consistency criteria 
(19) become independent of the pj-s: 

M MM 

p^Y^Pi^hM and \ = Y, Pik = hY^k . (27) 

i=l i=l i=l 



Substituting for p in the definition a — — (£ + jz^) and setting it equal to zero 
we obtain 

(1 - \)e = -hM. (28) 
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Fig. 2. Plot in (e, //)-space of the a = contour for M = 5 and li = i. 



Using the definition b = e^(l — A) we arrive at 

a" = ~ (29) 



We can see that for a given value of e < there is always a value of /i 
and consequently c below which the proteins spread out and the system is 
characterized by a low surface density but high surface coverage (a > 0), and 
above which the surface density is high but the proteins are relatively loosely 
adhered (a < 0) (see Fig. 2 as an example). 

Considering the number and size of the subsequent footprints immutable, 
the equilibrium state of the adsorbed proteins can be completely determined 
for any given value of e and fi. We can then plot the values of a, p or X in 
the [e, /i) parameter space and visualize the effect of varying the chemical po- 
tential (i.e. changing the concentration) or the adhesion energy (i.e. changing 
either the adsorbing protein or the surface [16,17]) for a given systemO] The 
change of sign of a in the (e, /i)-space corresponds, as shown in Eq. (29) to the 
curve /i = ln(— -^) (Fig. 2). This result is independent of the choice of {/j}-s. 
An examination of the unit surface coverage A shows a non-monotonic fi (or 
In c) vs. A dependence (Fig. 3 b, d, and e) characterized by a maximum and a 
minimum, while the total surface density p demonstrates similar unexpected 
non-monotonic behavior exhibiting a maximum in e (Fig. 3 a, c, and e). Ex- 
perimental evidences for this latter behavior can be found in Ref. [18], where 
the saturation level of the adsorbed amount of various proteins (HSA, hFb) 
showed a non-monotonic dependence on the grafting ratio of PEG chains on 
Ti02 surfaces. 



^ The chemical potential as defined by Eq. (15) also depends through /xq on e as 
well as fed and fea which probably depend on e as interpreted above, but a pure In c 
vs. e plot can be recovered from the p-o + lnc vs. e results by a transformation of 
the coordinates. 
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(a) (b) 




Fig. 3. (a) and (b) contours of A and p in the (e, /i)-space, (c) and (d) slices of A 
along /i at e = —1.2 and p along e at = 1, (e) surface plot of A and p together, 
for M = 5 and U = i throughout. 

2.4 Non-linear potentials 

For potentials with non-linear footprint size dependence, the equilibrium state 
density distribution can be calculated numerically just as easily as for the 
preceding linear case. The solution, on the other hand, is no longer a simple 
exponential distribution and the potential cannot be as easily parameterized 
in a physically meaningful manner. 
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3 Dynamics 



In the following we introduce two limiting cases of a general dynamics, both 
of which minimize the free energy density (16). In these models wc consider 
thcrmodynamically consistent rates for the transitions of a protein between 
neighboring states and for the adsorption and desorption of the protein from 
and to the smallest footprint state. These rates are derived from a potential 
for which each state is separated by a potential barrier (Fig. 1). Consequently, 
the rate of a protein expanding from state i (for 1 < i < M) to state i + 1 
(presuming there is sufficient adjacent empty space) is 

k+ = Uie 2fcBT , (30) 



while the rate of shrinking from state i (for 1 < i < M) to i — 1 is 

k- = Vi-i e 2feBT ^ (31) 



where Ai^j+i ^ = — and is determined by the height of the energy 
barrier. These rates satisfy detailed balance by construction. The rate of des- 
orption from state 1 is given as 

A;r = A;d, (32) 

while the adsorption rate per unit length (into state 1) is 

— kg^c. (33) 



To complete a dynamical model we must also have information on the available 
surface surrounding each protein at any given time as well as on the diffusion 
of proteins along the surface. 



3.1 The infinite diffusion limit 



If the diffusion of adsorbed proteins on the surface is considered infinitely fast 
(ideal mixing) we can construct a mean field-like model where every protein 
on the surface is presumed to "feel" the same environment. Defining the state 
density currents (i.e. the number of transitions per unit time and length) jf 
according to Fig. 4, we can see that every j^^ (corresponding to shrinking 
and desorption) is independent of the protein environment, whereas every jf 
(corresponding to spreading and adsorption) depends on the adjacent free 
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Fig. 4. The complete set of the state density currents . 

space, which is considered to be identical for each protein in our mean field- 
like model. These currents then have to include the probability of finding a 
sufficiently large empty interval to either side. Since we assume the proteins 
to be perfectly mixed, the distribution of empty intervals is exponential with 
a mean value of (L — l)/N = (1 — A)/p, whence the probability of finding an 
interval with a size of at least x is then 



P{x) 



e "^1-^ 



(34) 



Thus, we can write the state density currents as 



Ji+l — Pi+lK+1 — Pi+lK+1 

jt = PiktP{li+i - k) = PiH e-('^+-'^)T^ 



> i > 1, 



(35) 



and 



Ji = Pi fed = Pi fed 

Jo = PoKP{h) = Po/^ae"'^^ 



0. 



(36) 



Having derived the currents jf, our system can be completely described by 
the following M differential equations: 



Pi = -Ut + k) + {jti + k+i) i<i<M, 



(37) 



which can be readily solved numerically (as an example see Fig. 5a). 

To verify whether this dynamics converges to the equilibrium state distribu- 
tion, characterized by the minimum of the conditional free energy (16), let us 
calculate the rate of the change of this free energy if proteins in state i start 
to expand to state i + 1 (for 1 < i < M) : 



00 d(p pi+i 

^ T— = {Si+i - Ei) + In 

dpi+i dpi Pi 



+ 



P 



1-A 



{h+i — k)- 



(38) 
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This quantity is identical to 



In 



Ji+l 




+ ln 



Pi+i 



Pi 



+ 



1- A 



P 



(39) 



which means that expansion (j^) is faster than shrinking (jj+i) only if the rate 
of the free energy change for expansion is negative. A very similar derivation 
holds for the relation between the adsorption (jq") and desorption (jf) cur- 
rents. Thus indeed, the currents between any two states point to a descending 
direction on the free energy landscape, driving the system to equilibrium. 

3.2 The zero diffusion limit 

If on the other hand we assume the proteins adsorbed on the surface immo- 
bii^H after adsorption, we can no longer construct a mean field model such as 
above. Transient effects similar to jamming in RSA begin to play a significant 
role, the system relaxes very slowly, and the transient state distributions can 
be significantly different from the eventual equilibrium. To study the zero dif- 
fusion behavior of our model we conducted Monte Carlo (MC) simulations of 
systems of approximately 100 proteins (system size was fixed at lOOmaxj/,}) 
several hundred times and averaged the results. Slow relaxation and large 
differences in state distribution under similar total surface densities were ob- 
served (see Fig. 5) over a wide range of parameters. The simulations qualita- 
tively reproduced several known experimental effects, such as the overshot of 
the surface density at high concentrations or hysteresis of the surface density 
as a function of the concentration even after large equilibration times. We are 
performing experiments which will be analyzed by such MC simulations. 

3. 3 Exchange 

A further consequence of reversibility, aside from the subtle evolution of the 
state density, is the possibility of the exchange of the adsorbed proteins by 
the proteins of the solution (a phenomenon widely observed among simple 
homo-polymers [14]). To examine exchange we conducted simulations of a hy- 
pothetical setup were "labeled" proteins were allowed to adsorb for a specified 
time (2/z/) and then the solution was replaced by either a buffer solution (0 

^ Limited movement occurs none the less, through "crawling", i.e. successive 
spreading and receding of proteins resulting in movement of the particles' center 
of mass. 
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(a) 

total 




Fig. 5. (a) MC simulation of a protein-surface system with M = 3 possible states 
of sizes li = 1, I2 = 1.5, and I3 = 2 with binding energies £1 = —11.25 A;bT, 
£2 = —14.0625 ksT, and £2 = — 15 ksT, respectively. The states are separated by 
barriers of equal height, characterized by a state independent Vi = v prefactor in the 
corresponding transition rates, as defined in Eqs. (30) and (31). The desorption rate 
is also given by Eq. (31) with the hypothetical £0 assumed to be 0. The adsorption 
rate is set to k^ = vc with a dimensionless protein concentration of c = 0.2. The 
simulation was conducted with a system of size L = 200 and averaged over 10000 
runs. The state densities from the MC simulations (thick lines) as a function of time 
(measured in units of 1/z/) are plotted together with the equilibrium distribution 
(using Eq. (18)) (dashed lines) and the numerical solution in the infinite diffusion 
limit (Eq. (37)) (thin lines). The gray levels (starting from light gray) correspond to 
Pi, P2, P3, and p. (b) Sketches of different regimes of the dynamics: under low surface 
coverage all adsorbed proteins spread out to their energetically most favorable state 
(top row), while for high coverage - assuming a convex potential - the equally, but 
not fully spread out state may be preferential (bottom row). 

c unlabeled proteins), or an "unlabeled" solution of the same concentration 
(1 c unlabeled proteins), or a more concentrated unlabeled solution (1000 c 
unlabeled proteins). Observing the concentration of labeled proteins on the 
surface after the introduction of the new solution (Fig. 6), we can see that the 
density of labeled proteins decreases faster as the concentration of unlabeled 
proteins is increased. 

Further examining the various state densities, we see that this increase can be 
attributed to the labeled proteins being "ratcheted" up to smaller footprint 
states and eventually expelled from the surface by the newly adsorbed unla- 
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(a) 




1 2 3 4 5 



Time [ 1 /V ] 




Fig. 6. (a) Plot of the state densities of "labeled" proteins from the MC simulations 
of the system described in Fig. 5 (but averaged over only 2000 runs), using the same 
color code. As detailed in the text, the initially solution of labeled proteins (with 
concentration c = 0.2) was replaced at time 2/z/ by either a buffer solution (solid 
lines), or an "unlabeled" solution of the same concentration (densely dashed lines), 
or a more concentrated unlabeled solution (sparsely dashed lines). For comparison, 
we have also plotted the total surface density p in the infinite diffusion limit with 
buffer replacement (upper dotted line) , as well as in the limit of the fastest exchange 
(i.e., when no spreading of the labeled proteins is allowed) (lower dotted line), (b) 
Caricature of the "ratcheting" mechanism of protein exchange. 

bclcd proteins [15]. This is most obvious in case of the replacement with the 
most concentrated unlabeled solution (1000 c), where a large maximum of the 
surface density of state 1 (pi) appears near time l.hjv. Since desorption occurs 
exclusively from this state, the inflection point (i.e., the maximal slope) of the 
total surface density (p) near the same time is another manifestation of the 
same ratcheting phenomenon. Because this inflection at high concentrations is 
a clear sign of exchange and could easily be measured, we are currently setting 
up similar experiments with fluorescently labeled proteins. 



4 Discussion 

We have presented a thermodynamically consistent model of protein adsorp- 
tion possessing an equilibrium state. The existence of which facilitates the 



15 



understanding and prediction of long-term effects and also offers novel pre- 
dictions for the dependence of the equilibrium surface density (p) and surface 
coverage (A) on the binding conditions and the protein concentration. We have 
also predicted the existence of a well defined transition between predominantly 
compact and expanded states in the adsorbed layer as a function of either the 
concentration or the binding conditions. 

Introducing dynamics, we derived an analytically tractable mean field model, 
which is equivalent to the limit of infinite surface diffusion of the adsorbed 
proteins. We have also performed direct MC simulations of the zero diffusion 
dynamics, and have been able to qualitatively reproduce several experimental 
phenomena. The comparison between the mean field model and the MC simu- 
lations underpins the importance of RSA like surface frustration in the slowing 
down of the adsorption dynamics. The MC simulations have also led us to the 
idea of a novel experimental method for the study of protein exchange. 



A Appendix: Configurations on a continuous line 



Starting from a continuous model where each adsorbed protein is located along 
a continuous line (ID surface) of adsorption, approximation (8) can be shown 
to be exact in the (5 — >■ (no — > oo) limit. To arrive at this result, directly the 
continuous configuration volume C{{ni}) has to be derived. 

Let us first consider only one adsorbed protein (A^ = 1) in state i. This protein 
separates the adsorption surface into two disjunct intervalsPI each of which 
can have a length of maximum L — Zj. For arbitrary A^, each of the A^ -|- 1 
intervals has a length 

Wje[0,L-l], where j G {1, 2, .., A^ + 1} (A.l) 



with I being the total occupied surface. The intervals wj must also satisfy the 
spatial constraint 

Af+l 

J2wj = L-l. (A.2) 



This means that the configuration volume C{{ni}) for a single protein is simply 
the length of the hypotenuse of a right triangle with legs of equal length, L — I 
(Fig. A.la). 

^ This separation of intervals (or shielding property) only holds in ID, which is 
one of the main reasons that 2D calculation are far more difficult. 
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Fig. A.l. C({nJ) for (a) N = 1 and (b) iV = 2 

For N = 2, C{{ni}) is the area of the equilateral triangle base of the pyramid 
formed by the origin, (L — /, 0, 0), (0, L — /, 0), and (0, 0, L — /) in the 3D space 
of wi, W2 and ws (Fig. A. lb). So far we have considered the adsorbed proteins 
indistinguishable - by only considering the lengths and arrangement of the 
empty intervals, but not the order of the adsorbed proteins -, consequently 
this area has to be multiplied by a factor of two when the proteins are in 
different states. 

In general for N > 2 the configuration volume is an A^-dimensional hyper- 
surface in the 1-dimensional space of wj-s defined by the constraint (A. 2), 
multiplied by a factor of z=Jf^ — . This can be written in an integral form as 



°° N+l N+l N+l 



c{{n,}) = / n n ©K) s{l-i-y: (a.3) 



-oo 



where B(x) is the Heaviside function. Carrying out the integration with respect 
to wn+1 we have 



N+l N 

J dwN+i e{w^+^) 5(1-1- J2 ^j) = - ^ - E ^j)- (A-4) 

i=i i=i 



Integrating again with respect to wn '■ 



oo jy 

I dwN Qiw^) QiL-l-J2' 

-oo J=l 

N~l N-1 

= {L-l-J2w,) 0(1-1- 



(A.5) 



Now we can see that by integrating with respect to the remaining N — 1 
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variables one by one in a similar way, we arrive at 



CiM) = ^-^M-^, (A.6) 

which, divided by the unit volume , is the same as Eq. (7) with the approx- 
imation (8), which can now be considered exact in the 5 — > limit. 
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